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SYiTOPSIS 


The present vjorlc deals with m-raachine n-joh flow shop 
scheduling problem considering the total number of tardy 
jobs as the effectiveness criterion. Pour heuristic algorithms 
are proposed and have been tested, on 60 randomly generated 
problems of 12 different sizes. The performances of the 
heuristics have also been tested for a special class of 
ordered flow shop problem , ill the four heuristics gave the 
same results for the above special class of the ordered flow 
shop problem, considered in the present work. The performances 
of the heuristics were evaluated based on average number of tardy 
tardy jobs, number of times best solution was generated, the 
frequency of obtaining guaranteed optimal solution and the 
average execution time for various problem sizes. The analysis 
indicated that, heuristic-4 gives overall best performance. 
Further, it was observed that the number of jobs have more 
significant 'effect on the computational times than the number 
of machines. In case of all the heuristics the average exe- 
cution time is veiy small. It is suggested that a problem 
should be solved using all the four heuristics and the best 
solution should be adapted since no heuristic explicitly 
guaranties the best solution. 
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INTROBUCTIOIT 

Scheduling can he defined as the allocation of resources 
over time to perform a collection of tasks. lii a typical 
production situation, where the resources are generally 
machines and the tasks are jobs, scheduling decisions are 
generally reached by using the scientific and systematic 
approaches involving the following four steps, 

1 ) Pomulation 

2) Analysis 

3) Synthesis 

4) Evaluation 

In the first stage, basically a problem is identified 
and the criteria that should guide decision-making are 
determined. This is often a subtle and complicated activi- 
ty, but good decisions can seldom be expected without a 
clear definition of the problem at hand and an explicit 
recognition •* of objectives. Analysis is the detailed pro- 
cess of examining the elements of a problem and their inter- 
relationships. This stage is aimed at identifying the 
decision variables and also at specifying the relationships 
among them and the constraints they must obey. Synthesis is 
the process of building alternative solutions to the problem. 
Its role is to characterize the feasible options that are '■ 



2 


available, finally, evaluation is the process of comparing 
these feasible alternatives and selecting a desirable course 
of action. This selection, is, of course, based on the 
criteria that are developed afc the outset. 

Ideally, the objective function should consist of all 
the system costs that depend on scheduling decisions. In 
practice, however, such costs are difficult to measure, 
identify and incorporated completely. Moreover, objective 
functions consisting of more than single type of costs (which 
may be contradictoiy in nature) sometimes become too 
complex to be analysed. Broadly, three types of decision 
making goals seem to be prevalent in scheduling : (i) effi- 
cient utilization of resources, (ii) rapid response to demands 
and (iii) close confomance to prescribed dead lines, Pre- 
quently an important cost related measure of system perfor- 
mance (such as, machine idle time, job waiting time, job 
lateness, etc.) can be used as a substitute for the total 
systan cost, and quantitative approaches to problems with 
these criteria appear throughout the literature on scheduling. 

Two kinds of constraints are generally encountered in 
scheduling problems, First, there are limits on the capa- 
city of available resourbes, and, second there are technolo- 
gical restrictions on the order in which tasks can be per- 
formed . 

iny solution to scheduling problem has to answer the 
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tv/o major questions. 

1) Wliich resources will be allocated to perform each task? 

2) When will each task be performed? 

In other words, the essence of scheduling problem gives 
rise to allocation and sequencing decisions, Ihe present work 
deals with the multistage flow shop scheduling problem. In 
flow shop jobs are multistage in nature (a job is a set of 
operations with the special precedence relations) and machines 
have serial configuration, i.e,, it is possible to number the 
machines so that if the ith operation of any job j precedes 
its kth operation then the machine required by ith operation 
has a lower number then the machine required by the kth ope- 
ration. In particular each operation after the first has only 
exactly one direct predecessor and each operation before the 
last has exactly one successor. In general any flow shop 
scheduling problem can be defined as follows : 

Siven a flow shop environment with n jobs and m machines 
the objective is to determine an optimal sequence considering 
some measure of effectiveness. The effectiveness of any 
sequence can be measured in terms of the criterion of effe- 
ctiveness, Some of the important effectiveness criteria 
are : make-span time, average completion time, due date per- 
formance, machine utilization, inventory of jobs in process 
etc, Ihe due date performance of any sequence can be measured 
in terms of lateness of the job, total or mean tardiness. 
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maxiFiiiKi "tardiness and ■to'tal number of "tardy jobs. Iia"teness 
of a job can be defined as "the difference be"t"ween i"ts com- 
pletion time and its due date, The maximum of labeness and 
zero is called tardiness. A job is called tardy or late if 
its tardiness is non -aero (i.e., strictly positive) and non- 
tardy or early otherv/ise. 

Considerable effort has been directed by many researchers 
towards the solution of flow shep scheduling problem, subject 
to standard set of assumptions listed by Bakshi et al [1]. 

Most of the flow shop problems have been tackled using com- 
binatorial optimization techniques, Monte Ca.rlo sampling, 
random sampling, biased sampling, decomposition, heuristic 
approaches and complete enumeration techniques. Most commonly 
used combinatorial method is the Branch and Bound methodology. 
Minimization of the make span time and average flow time are 
the criteria considered most often. Jhonson 12] has suggested 
an algorithm which yields optimal solution for two machine 
n job flow shop case considering minimization of the make- 
span time. Number of lieuristic approaches, have been reported 
in the literature, for the general mxn flow shop scheduling 
problem considering minimization of make span time as the 
effectiveness criterion. Some of the notable contributions 
are due to Campbell Budek aiid Smith [ 31 , Palmer [4 J , G-upta, 
J.N.D {51 , Petrov {61, Nobesimha {71, Krone et al {8}, Page (9 } 
etc. Other criteria like minimizing the total or mean tardiness 
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minimizing the inventory of johs in process, maximizing the 
machine utilization etc. have also hc-en considered. Survey 
of the literature suggested that vexy little work has keen 
done for general flow shoo scheduling problems considering 
due date related criteria. In particular no algorithm has 
been reported for the minimization of the total number of 
tardy jobs. May be, this is due to the fact that the per- 
fornic^ice of algorithms considering the due date related cri- 
teria is not linearly related to the completion time as the 
value of the objective function depends upon two independent 
entities, viz., the completion time of the job (which is 
senuence dependent) and Its due date (which is independent 
of the job sequence). Hence, in general, the analysis of 
the due date related criteria is more difficult as compared 
to other criteria. 

Maxwell 1101 has given an integer programming formula- 
tion for one machine n job flov; shop problem to minimize the 
total number of tardy jobs. This formulation was questioned 
by Sidney C11l and he gave a corinter example to show that the ; 
cuts (cutting planes) generated by Max^\'ell's procedure elimi- 
nate the optimal solution. Maxwell f12j has acknowledged | 

that this algorithm does not generate the efficient cuts. For | 

i 

the counter example given by Sidney, he has shown that his i 

procedure leads to an alternate optimal solution. Further, | 

I 

[ 

f 

! 

I 

i 

i 
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Moore [13 1 claims that the proof given hy Maxwell is incom- 
plete and the anraher of variable? and constraints expand 
rapidly and hence for the case of multistage flov; shop 
problem this approach seems infeasible. 

Moore [ 14] h.as given a simpler met hod for 1xn flo'w shop 
problem for the criterion of minimization of total number of 
tardy oobs. Many investigators [131, [161, fill, [14* have 
given proof of optimality of Moore's algorithm or its closely 
related form. 

In the present -work four heuristic approaches have been 
developed ior the general mxn flow shop scheduling problem 
considering the criterion of minimization of the total number 
of tardy jobs. The conceptual frame work suggested by Moore 
for single machine case has been utilized effectively for 
the development of these heuristics. The proposed heuristics 
have also been tested for a special class of mxn flow shop 
scheduling problems, known as the problem with ordered 
processing time matrices. 

The thesis is organised in three chapters. In Chapter 
II problem statement and solution methodology are presented. 
Computational experience with large number of randomly 
generated problems of machine size varying from 2 to 5 and 
the job size varying from 10 to 50 is presented in Chapter III 
alongwith the conclusions. 
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CHiPTER 2 

PROBLEM PORMULiPIOF MD SOLULIOE 
METHOLOLOGY 


2.1 SL^EMEET OP THE PROBIEM 

Consider a multi-machine flov shop situation with m machines 
and n jobs. Each job is made up of a set of tasks or operations 
which are to be processed in a given technological order, ill 
the standard assimiptions given in [ 1 ] for a general flow shop 
environment are considered as valid and the objective is to mini- 
mize the total number of tardy jobs (H^). 

2.2 HOMEHOLiTURE 


m 

n 

( j) 

■t(i, j) 
d(j) 
c(i, j) 

% 

- {E } 
{I } 
L(j) 
T(j) 


Humber of machines 

Humber of jobs 

jth position in any sequence 

- Processing time of job (j) on machine i 
Due date of job (j) 

- Completion time of job (j) on machine i 

Total number of tardy jobs for the given sequence 

Set; of early or non-tardy jobs 

Set of late or tardy jobs 

Lateness of job(j); L(j) = c(ra,j) - d(j) 

- Tardiness or positive lateness of job (j); 

T(j) = max {0, L(j)} 
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IT(3) - Difference of the completion times of the joh (j) 

and job (o-1) on the last machine; IT (o) = 0(1,3) 
C(i,3-1) 

S(d) - Sum of the processing times of (j) on all m 

machines; 

m 

S(3) = I t(i,3) 
i=1 

2.3 MlTHMlTIGili DOmULATIOF 

Mathematically, the problem can be stated as follows ! 
n 

Minimize = I P(Tjp(3)) 

0=1 

where DCx) = 1 , if x > 0 
= 0 otherwise. 

Subject to the set of standard flow shop constraints as listed 
by Bakshi et al [11, 

2.4 FOUND jffilOFS OF PROPOSED METHODOLOOr 

For the case of single machine, n job flow shop scheduling 
problem, Moore [ 141 proposed an algorithm for the minimization 
of the total number of tardy jobs. The algorithm alongwith 
the proof is s'fltemarized below. Moore' s algorithm is based on 
the following result, the proof of which is self evident. 

In any situation, if the EDD (early due date) sequence 
yields zero or exactly one tardy job, then it is the optimal 
sequence for the number of tardy jobs criterion. However, if 
it yields more than one tardy job the EDD sequence may not be 
optimal. 
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The algorithm suggested hy Moore [ 14 1 assumes the form 
of the optimal sequence as follows ; 

a) First, a set {E* } of early or non-tardy ^obs in EDE order 

b) Then, a set {1* } of late or tardy jobs in any order* 



Figure 1 The form of a sequence that minimizes 

2,5 MOORE'S iLGORITHM 

The various steps of the algorithm are : 

Step-1 : Place all the jobs in set E according to EKD order. 

Let set L be empty. 

Step- 2 ; If no jobs in set E are late, stop: E must be optimal, 
otherwise identify the first late job in set E. 
Suppose this turns out to be job (j). 

Step-3 : Identify the longest job among the first j jobs in 
sequence. Remove this job from set E and place it 
in set L. Revise the completion times of the jobs 
remaining in set E and return to step-2, 



2,5.1 Proof : 


As the job (d- 1) is non-tardy, its due date is greater 
than or equal to its completion time and as the 30b (3) is 
tardy, its due date is less than its completion time. Hence, 



0(3-1) < 

<i(j-1) 

(1) 

and 

0(3) > 

<i(j) 

(2) 

but 

0(3-1) = 

3-1 

y t(i) 

(3) 



i=1 


and 


i=1 



0(3) = 

I t(i) 

(4) 


i=1 

Substituting (3) and (4) in (1) and (2), -we get, 


r t(i) y d(3) 
i=1 

or 


3-1 

^ t(i) + t(3) > d3 (5) 

i=1 


and 


3-1 

I t(i).^ d(3-1) (6) 

i=1 


After the job with the maximum processing time among the 
first 3 jobs is moved to set 1, the completion time of job 
(j) would be reduced atleast by the amount t(j). Otherwise, 
if max t(i) occurs for i = j, then job (3) would be moved to 

1 < i < j 

set I. In the latter case the number of tardy jobs does not 
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increase as the joh (j) -was tardy anyway. In the former case, 
also the number of tardy jobs does not increase. Since in 
this case the completion time of job (j) is reduced by 

t(k) = max {t(i) } As t(k) is greater than t(j), its new 

completion time becomes less than the left hand side of eqn. 

(6), which is less than or equal to d(j-1) and hence job (j) 
becomes non-taidy, 

2.6 PROPOSED MElHODOlOOr 

2,6.1 Development 

As stated earlier in Chapter 1 the IP approach of Maxwell 
{10] has two main drawbacks even for 1xn problem, firstly , the 
cuts (cutting planes) generated are very inefficient and there- 
fore the greater computational effort is required. Secondly , 
the number of variables and constraints beoomesprohibitively 
large as n increases. Thus, as such the IP approach has hardly 
any use even for the moderate size real life mxn flow shop 
problems. Hence, the only feasible alternative left is to go for 
a heuristic approach procedure which may yield fairly good solu 
tions without requiring large computational effort, 

Moore's algorithm I 141 for Ixn problem is fairly simple 
and does not require large computational effort. As pointed 
out in Section 2.5, the algorithm requires the lenghiest job to 
be identified. The purpose of identifying the lengthiest job 
is that it would be removed from B and would be placed in set D 
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in order to reduce the completion time of the first tardy job 
to an extent that it becomes non-tardy, for a single machine 
case, it is assured that by shifting the lengthiest job to set 
L, atleast the first tardy job -would become non-tardy as the 
lengthiest job is uniquely defined due to the absence of idle 
time, for a multimachine case, all the machines but for the 
first machine -will generally have inserted idle times. The 
completion time of any job is made up not only of the processing 
time, but the idle time is also included in it. Hence, the job 
to be shifted to set 1 from among first j jobs -where job (j) 
is the first tardy job should be the one which brings in maximum 
reduction in the completion time of the first tardy job. let 
us term this job as the 'shift' job. But the main difficulty 
here is that the idle times are sequence dependent and would be 
changed if the sequence is changed. As such^one can not guarantee 
that by shifting any particular job to set 1 the first tardy job ; 
would become non-tardy. Therefore, the strategy to choose the : 
shift job to be shifted to set 1 is very impoirtant. The value 
of the objective function is calculated after this job is shifted| 

f 

to set 1 and compared with the objective function value before 

i 

the change in the sequence was made. The sequence giving the ; 
better value of the objective is retained. ! 

f 

Thus, the important differences betweaa Moore's I 141 singl| 

machine algorithm and the proposed heuristics are : (1) the value] 

i 

of the objective function is to be calculated in the proposed | 



methodology "before and after the shift job is shifted to set 
L and only if the proposed change in the sequence brings in 
some improvement in the objective function value, it is to be 
carried out. Another important difference is the vioy in -which 
the job to be shifted to set L is chosen. Pour heuristic- rules 
are developed for selecting the shift job, -which form the basts 
of the proposed heuristics. 

2,6.2 Heuristic ~ 1 ; 

The principle aim in selecting the shift job is to see that 
the removal of the shift job from E brings in the maximum attai- 
nable reduction in the completion time of the first tardy job 
and makes it non-tardy. Suppose job (j) is the first tardy job, 
let 

m 

S(k) = J t(i,k) denote the sum of ths processing tii 
i=1 

on all m machines for any job (k). Select the shift job as the 
one having the maximum value of S(k) for k = 1,2, j* It is 

hoped that by choosing the job having maximum S(k) as the shift : 
job and removing it from B would make : 1) the job (j) non-tardy 
and 2) its new completion time not greater than the previous ; 
completion time of job (j-1). However, it is not assured that 

I 

the shifting operation i.e., removal of the shift job from E 
would always accomplish objectives (1) and (2). Therefore, i 

Hm is calculated before and after the shifting operation. If i 
the proposed shifting operation increases job (j) becomes I 
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the new shift job instead of job (k) and the shifting operation 
is carried out. If, however, objective (1) is accomplished 
and Fij, remains the same by the proposed shifting, then the 
accomplishment of objective (2) is checked. If objective (2) is 
also accomplished then the proposed shifting operation is per- 
formed! otherwise job (j) becomes the shift job instead of job 
(k) and the shifting operation is carried out. However, if 
objective (1) is accomplished and Hg, decreases by the proposed 
shifting, job (k) stays as the shift job and the shifting opera- 
tion is carried out. 

The methodology mentioned above for selecting the shift 
job is the basis of the present heuristic, ill other steps are 
similar to Moore's [ 14 ] single machine algorithm. 

The various steps of the algorithm are : 

Step-1 ; Place all the jobs in E according to EDD order. 

Let set L be empty. 

Step-2 ! If no jobs in E are late, stop; otherwise, identify 
the first tardy job in E. Let job (j) be the first 
tardy job. Go to step-5. 

Step-3 : Identify the job having the maximum value of S(i) 
for i = 1,2, ... j. Suppose this is job (k) 

Go to step 4. 

Step-4 : If j = k, remove this job from B, place it in L and 
go to step-5. Othersie, go to step-6. 
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Step-5 : Revise the completion times of all the jobs in S. 

G-o to step-2, 

Step-6 : Calculate Call it . Calculate Rg, hy removing 

joh (k) from E and placing it in E. Call it R/ji2* 

There are three cases to be considered. 

Case (i) ; > F ^2 

Remove 3ob (k) from E and place it in E and go to step-5. 

Case (ii): Fg,-] <- F^ 2 ' 

Remove job (j) from E and place it in E and go to step-5. 

Case (iii): Fg,^ = Fgi 2 

Calculate C(m,3-1). Call it C(m,j-1)^, Calculate C(m, j-1) 
by removing job (k) from E and placing it in E. Call it C(m,j-1)2* 
If 0(m,j-1)^ 2 C!(m,j-1)2, remove job (k) from E and place it in 
L, Go to step-5. Otherwise remove job (j) from E and place it in 
L. Go to step-5. 

2,6.3 Heuristic- 2 ; : 

Let us examine the general expression for the completion I 

time of any job (i), which is given below : | 

C(m,i) =r t(m,i) + max {G(m, i-1 ) , C(m-1, i)) i 

i 

From the above expression it can be clearly seen that the | 

term t(m,i) would always contribute to the completion time of | 

job (i). larger value of t(m,i) would result in larger C(m,i) | 
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and sulDsequent delay in the following jobs, therefore, if the 
job (k) having the maximum value of.t(m,i) for i = 1,2, ... j 
is selected as the shift job, it can be hoped that the shifting 
operation would accomplish objectives (1) and (2) stated in 
Section 2,6,2. As pointed out earlier, the proposed shifting 
is carried out only if it accomplishes objectives (1) and (2). 
Otherwise job (j), the first tardy job becomes the shift job in- 
stead of job (k) and the shifting operation is performed. This 
methodology of selecting the shift job is the basis of the 
present heuristic. The various steps of this algorithm are identic! 
that of to ■ heuristic-1 except step-3. Step-3 of the heuri- 
stic algorithm 2 is given below: 

Step 3 : Identify the job, having the maximum value of t(m,i) 
for i = 1,2, ...'j. Suppose this is job (k). Go to 
step 4. 

2,6,4 Heuristic - 3 : 

Por any sequence the processing times of the jobs on the 
first machine play an important role to determine the completion , 

times of the jobs. If a job has a larger processing time on 
the first machine, the subsequent machines have to be idle, ' 

waiting for the job to arrive from the first machine. This i 

waiting time, in turn increases the job completion time and all ; 
the subsequent jobs are delved. Hence, in order to reduce the S 
waiting times and in turn completion times of all the subsequent | 
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jobs in the sequence, the shift jot should be the one having 
the maximum value of t(1,i), for i = 1,2, j, where job (j) 

is the first tardy job. let t(1,k) be the maximum among first j 
jobs. Then with job (k) as shift job, the accomplishment of 
objectives (1) and (2) stated in Section 2.6.2 is checked. If 
both are accomplished the shifting is performed. Otherwise, job 
(j) becomes the shift job and the shifting operation is performed. 
This rule of selecting the shift job forms the basis of the 
present heuristic. The various steps of this algorithm are 
identical to heuristic - 1 except step 3. Step 5 of the heuri- 
stic algorithm 3 is given belov/ : 

Step-3 : Identify the job having the maximum value of t(1,i) 
for i = 1 , 2, ... j. let this job be (k), G-o to 
step- 4. 

2.6.5 Heuristic - 4 : 

The information of the idle time for the last machine, i.e., 
time span during which the last machine is idle and waiting for 
a job to be processed is very important. The idle time on the 
last machine for any job reflects the effect of its processing 
times on all the m machines on its completion time. If a job 
has a larger processing time on anyone of the m machine it 
would result in a larger idle time on the last machine. Thus, 
the vital information of idle time on the last machine. "Ibe , 
point to] be borne in mind while exploiting this important and 
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useful information is that the idle time on the last machine 
does not reflect the effect of the processing time on the last 
machine. For example, there could he a joh in the sequence for 
■which the processing times on all the machines hut for the last 
machine are smaller in magnitude and as a result it produces 
very little idle time on the last machine. But, if, this job 
has a larger processing time on the last machine, the completion 
time of the job ■will he large and this "will never he reflected 
by the idle time on the last machine. Therefore, sum of the 
idle time and the processing time on the last machine of the 
various jobs is taken as the criterion to identify the shift 
job. In this heuristic, the 30b having the maximum value of the 
sum of idle time and the processing time on the last machine is 
selected as the shift 30b. The sum of the idle time and pro- 
cessing time on the last machine could be vie-wed as the modified 
processing time on the last machine, meaning that 30b (k) virtua- 
lly requires time equal to the sum of the above two quantities on 
the last machine. Either the last machine is waiting for the 
arrival of job (k) or processing it. If the above mentioned sum 
for job (j) is denoted by IT (j), then the expression to calcu- 
late it can be derived as follows ; Following figure gives a Gantt 
chairt for a two machine, four job problem. 
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Machine 1 


Machine 2 




'p;~^ iTTr T - r- rr: 


1 ( 2 ) 1 ( 3 ) 

IT(|kI0?(3) 

f " , I 

f 

Time 


IT (4) 


Pig, 2 Gantt chart for a 2-maohine 4 joh problem 


Prom the figure 2 given above, it can be seen that 

IT(2) = C(2,2) - C(2,1) = t(2,2) + 1(2) 

IT(3) = G(2,3) - 0(2,2) = t(2,3) + 1(3) 

IT(4) = 0(2,4) - 0(2,3) = t(2,4). 


Prom above vie can write a general expression for IT (2) as 
follows : 

IT(k) = 0(m,k) - 0 (m, 3 - 1 ) for k = 2,3, ... n. 

It should be noted that IT(1) is not defined because 0(iii,0) 
does Qot exist. Therefore, if the first late 30 b occurs for 

3 = 2 , then the accomplishment of objectives ( 1 ) and ( 2 ) is ; 

[ 

checked with job (1) as the shift job. If both the objectives 
are accomplished, the shifting operation is performed. Other- ; 
wise, job ( 2 ) becomes the shift job and the shifting operation j 
is carried out. If the first late job does not occur for j = 2 , ^ 
then the job having the maximum value of IT(k) for k = 2,3, ... j 
is selected as the shift job and then a check is made to see if ; 
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the objectives (1) and (2) are achieved. If both the objective 8 
are achieved job (k) stays as the shift job, otherwise job (j) 
is selected as the shift job in place of job (k) and the shifting 
operation is carried out. The procedure mentioned, above of 
selecting the shift job forms the basis of the present heuristic. 
The various steps of this algorithm are identical to heuristic-1 
except step 3, Step 3 of the heuristic algorithm-4 is given 
below t 

Step-3 ; If j = 2, let k = 1; go to step-4 j otherwise if 

j > 2, then identify the job having the maximum value 

of IT (k), for k = 2,3, ... j. Let this be job (k). 

Go to step-4. 

The common flow chart for all the heuristic algorithms is given 
in Appendix A. 

2.8 ILLUSTR1T3D EXiMPLl 

The suggested heuristics are illustrated with the help of 
a 2 machine, 5 job problem. The data for the problem is given 
in Table-I. Step-wise functioning of the various algorithm is 
given below : 

2.8.1 Solution by Heuristic-1 

The step-wise functioning of the algorithm is given in 


Table - II. 



Table II : Step-wise progress of Heuristic-I 



STAGS- 1 


step-1 

Initialize. E = {1-2-3-5_4} ; L = 0 


St ep-2 

Job 3 is the first tardy job in E 


St ep-3 

S(3) has the maximum value among the jobs 

in E 

St ep~4 

Job 3 is removed from E. S = {1-2-5-4 ) ; 

E = { 3 } 

Step -5 

Completion times of all the jobs in S are 

revised 


STAGE- 2 


St ep-2 

Job 4 is the first tardy job 


St ep-3 

S(4) has the maximum value among the jobs 

in E 

St ep-4 

Job 4 is removed from E. E = { 1-2-5 ) ; I 

= { 3,4 } 

Step-5 

Completion times of all the jobs in E are 

revised 


STAGE-3 


St ep-2 

ITo jobs in E are tardy. Solution is, 



{ 1-2-5-3-4 } or { 1-2-5-4-3 > and 

% = 2. 


In Stage-1, the jobs are arranged by EDL order and job 5 is 
foimd to be the first tardy job. The shift job upto and including 
job 3 is job 3. Therefore, it is shifted to L. Similarly, at 
Stage-2 job 4 is shifted to L. At Stage-3 no tardy jobs remain 
in E. Therefore, the algorithm yields two alternate sequences, 
viz. , 

1-.2-5-5-4 and 1-2-5-4-3, as the jobs in L can 

-flit 
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appear in any order. lor both the above sequences, = 2. 

2.7.2 Solution by Heuristic-2 : 

The various steps involved in the application of this 
heuristic are the sane as that of heuristic-1 except Step-3. 
Applying Stcp-3 of this heuristic it is found that the same 
jobs as sclccl ed by hcuristic-1 are chosen as shift jobs viz,, 
job 3 and 4 in stage 1 and 2 respectively. In St age- 3 F tardy 
jobs remain in E. Therefore the algorithm yields the same 
alternate sequences as that of heuristic-1. The step-wise 
progress of the algorithm is same as that of hcuristic-1* which 
is given in Table II. 

2.7.3 Solution by Heuristic-3 : 

The various stops involved in the application of this heuri- 
stic are the same as that hcuristic-1 except step-3. Applying 
step-3 of this algorithm it is found that job 2 has the maximum 
processing time on the first machine. It is removed from E and 

Fm^ , Fmo, C(m,j-1)-, and C(m.,j-1)p are calculated. As 0(m , 3 - 1 ) 2 ^ 

C(m,j-1)^ 

job 3 becomes the shift job and is shifted to set L. Stage-2 
and Stage-3 are identical to that of heuristic-1. 

The solution obtained by this heuristic is the same as that 
obtained by the other heuristics. Step-wise progress of this 
algorithm is given in Table U-A. 
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Processing times and due date data for the illustrated 

example 


■^1k 

■^2k 


Job k 

1 2 3 4 5 

1 2 4 9 7 

5 9 4 7 8 

5 13 14 20 22 


TABLE II-A 
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Step-wise Progress of Heuristic-3 


ST AGE- 1 

Step-1 and 2 are the same as that of Table II 
Step-3; t(1,2) has the maximum value among the job in E 
StGp-6; j%i2 3^*3. C(m,j-1)2 are calculated by removing job 2 from E 
■^Tl “ ^T2 D-'1)2* 3!herefore job 3 

is shifted to set L, E = {l -2-5-4) aid 1 = {3). 

3tep-5: Oompletion times of all the jobs in E are revised. 

Stage-2 and 3 arc identical to that of Table II. 

The solution is 1-2-5-4-3 or 1-2-5-3-4. 


2.7.4 Solution by Heuristic-4 

The various steps involved in the application of this 
heuristic are the same as that of heuristic-1. By applying 
step-3 of this algorithm it is found that, by chance, this 
heuristic also selects the same jobs as the shift jobs as 
selected by heuristic-1 at each stage. 

Job 3 and job 4 get selected as the shift jobs in stage-1 
and 2 respootivoly . At stage— 3 no tardy jobs remain in E. 
Therefore the algorithm yields the same solution as that of 
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hcuristic-l . The step-wise progress of the algorithm is same as 
that of heuristic-1 which is given in Table II. 

2.8 iPPIilCiTIOF OP VilRIOUS HEURISTICS TO THE ORDERED ili0¥ SHOP 
PROBIBH 

inothor class of scheduling problems v.;hich has attracted the 
attention of researchers recently is the mxn flow shop problem 
with ordered processing time matrices. This problem has two 
unique characteristics. 

i) If a particular job has a smaller processing time on any 
machine than any machine than any other job, then the 
procGssiing time of the former job will be less than or equal 
to the processing time of the later job on all corresponding 
machines. 

ii) The machine with the minimum processing time for a given 
job will also have the minimum processing time for every 
other job. 

Smith ot al [19 1 have suggested a methodology which gives 
the optimal solution considering the make-span as the effectiveness 
criterion for a special class of the above problem for which the 
maximum processing times occur either on the first or the last 
machine. 

In the present work, the heuristics developed in Section 2,6 
arc applied to a special class of the ordered flow shop problem 
for which the following additional constraints are valid. 
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) Proceosing times of any jot) k are in non-decreasing order 
of magnitude i.e., t(1,k) < t(2,k) $ 4 t(m,k) 

2) min { t(j,k) } .•> rjiax { t(i,k)} for i < j 
^ ^ ^ 1 .< k n 

As the processing times are in nondecreasing order, any 
machine will not have idle time. Therefore, if (j) is the first 
tardy job, then, 

C(i,j) = C(i,j-1) + t(m,j) > d^. ( 7 ) 

0(i,j-1) .< d(j-1) 

m-1 j-1 

hut C(i,j-1) = I l(i,1) + I t(m,k) . (8) 

i=1 k=1 


m 3-1 

I t(i,1) + I t(m,k) + t(m, 3 ) > dCj) (9) 

i=1 k=1 

Substituting for 0(i,3-1) in eOn. (^) and (g), we get, 
m 3-1 

I t(i,1) + r t(m,k) ^ d(j-1) (10) 

i=1 k=1 


The job having the maximum value of t(m,k) for k = 1,2, 

... j, is chosen as the shift job. It is assured that after 

the removal of the shift job from E atleast job (j) would 

become nontardy as its completion time would be reduced by 

max { t(m,k)}; sotherwise if the maximum occurs for i = j, 

1 ^ k < 3 

then jdb (j) would be shifted to L. Because of the special 
structure of the problem all the four heuristics would choose 
the sane shift job and hence proceed identically. Further tho 
special structure of the problem results in absence of idle times 
and thus all the four heuristics will yield optimum solutions. 
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CHiPTER 3 

RESULTS MD UISOUSSIOIS 

A composite computer code v;as developed for all the four 
heurirtios descriUed in Chapter 2. The proTramme v;as •written 
in I’ORlRAN’-lO and implemented on LEO-1 090 time sharinn computer 
systera. The computer listing of the programme is given in 
Appendix B. 

The performances of the various heuristics "were evaluated 
for 60 randomly generated prohlems of 12 different sizes. The 
number of machines considered were 2,3,4 and 5 while the number 
of 3obs considered were 10,25 and 40. Bor each problem size 
5 randomly generated problems were considered. 

5.1 GEiCSRAIIOI OB RMLOM PROBLEMS 

Bor each job, its due dates and the processing times on 
various machines were generated randomly from a uniform dis- 
tribution. The uniform distribution was assigned an upper 
limit on the processing times of jobs on various machines for 
each problem size. Similarly, the uniform distribution for 
determining the due dates was given an upper limit. The 
selected upper limits are given in Table 3. The lower limits 
for both due date and processing times were kept at zero. 



3.2 RESULTS 


3.2.1 PerfoiTianoe of Heuristics for Geaeral Plow Shop Problem 

Table 4 gives the relative perfomarce of the various 
hfurintiCB, It is observed that heuristics-1,2,3 and 4 yielded 
the bent solutions v.lth a freqi^ency of 7,14,1? and 17, res- 
pectively. This indicated that heuristic-3 and 4 are margi- 
nally ^'uperior to heuristic-2. Purther, the performance of 
the heuristic-l to generate the best sequences was inferior 
radongst all the four heuristics considered. Of the 60 pro- 
hlc no considered, all the four heuristics resulted in same 
setJuea^GS for 15 problems. 


Table 5 gives the relative performance of the various 
liouri-^tics and EI)I) considering number of problems for v/hich 
thi' optimality of the sequence is guaranteed (17^ = 0 or 1). 
Heurirt ics-1 and 2 resulted in guaranteed optimal solution for 
tlirec* additional problems while for ■hcuristics-3 and 4 the 
number was 4 additional problems. 

The relative performance of the various heuristics con- 
sidering the average number of tardy jobs (average basod on 
the five problems of the same size) is given in Table— 6. 

At the bottom of the table for each heuristic, the overall 
average number of tardy jobs considering all the 12 problem 
sizes is also given. It is observed that based on the 
criterion the performance of the heuristic 4 is the best 
followed by heuristic-2,3 and 1. 
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The larger size problem took about 102 milliseconds of 
execution time. 

Tt ’.var observed that the computational time increases 
r;iiinifiC''intly vjith an increase in the number of 3obs keeping 
the number of machines constant. However, the increase was 
not significant as the number of machines were increased keeping 
the nuinher of jobs 


constant. 
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3.2.2 Performance of Heuristics for Ordered Plow Shop Problem 

Tho performance of all the four heuristics was tested on 
3 orderer] flow shop problems given in Tables 10-12. These 
proolctfls carry the special structure, namely; 1) the processing 
time ol all the jobs on machine 1 is less than the processing 
tirae of tho same job on 2nd machine in the sequence; 2) the 
minimum processing time on the 2nd machine is greater than or 
equal to the maximum processing time on machine 1. 

Por each problem, all the four heuristics yielded the 
Simie results. The number of tardy jobs for problems 1,2 and 3 
were 1,4 and 4, respectively. The application of EDD rule gave 
2,4 and 6 as the number of tardy jobs for problems 1,2 and 3, 
respectively, 

3.3 CO INCLUSIONS 

In this thesis, four heuristics are developed for the 
general mxn flow shop scheduling problem for the miaimimization 
of the total number of tardy jobs. The performances of these 
heuristics were checked for randomly generated problems of 
12 different sizes. For each problem size, 5 problems were 
considered, The performances of the heuristics were evaluated 
based on average number of tardy jobs, number of times best 
solution was generated, the frequency of obtaining guaranteed 
optimal solution and the average execution time for various 
problem sizes. The analysis indicates that, heurist3.c-4 gxves 

overall best performance. 
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As 1 b cor et i cally envisaged, all the four heuristics 
yielded, same results for the 3 specially structured ordered, 
processing times problems. 

I'urther, it was ohseived that the number of jobs have 
more significant effect on the computational times than the 
number of machines. In case of all the heuristics the average 
execution ti.ue is very small. It is suggested that a problem 
should be solved using all the four heuristics and the best 
solution should be adapted since no heuristic explicitly 
guarantees the best solution. 
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TilDle 3 


Ilpi'cr 

Units on Prooeooin? Times and 

Due Pates 



r Various Prolnlcm i^izcs Oonsidored 


Pro])lum 

Uppp-r limit 

Upper limit 


ei7»c 

on piT.cerning 

on due dates 



tines 



?xlO 

30 

250 


2x23 

10 

250 


2x40 

20 

700 


3x10 

30 

450 


3x25 

10 

300 


3TC40 

20 

900 


4x10 

30 

750 


4x25 

10 

450 


4x40 

20 

1600 


5x10 

30 

750 


5x25 

10 

600 


5x40 

20 

1300 
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Table 4 


Relative Perforais’nce of Heuriatics 


Problem No. of 
sir.c problems 

mxn conside- 

red 


No. of problems for wbicli 
the best solution was 
obtained by heuristics No. 

12 3 4 


ITo. of problems 
for which all the 
heuristics gave 
ihe same solution 


2x10 5 20 

2x25 5 0 1 

2x40 5 1 5 


2 

2 

1 


3x10 5 

3x25 5 

3x40 5 


1 1 0 1 
0 1 1 2 
1 0 2 0 


1 

3 

2 


4x10 5 

4x25 5 

4x40 5 


0 1 
1 0 
0 2 


1 1 

3 1 

1 1 


1 

1 


5x10 5 

5x25 5 

5x40 5 


0 1 
0 2 
1 2 


0 1 

2 3 

5 1 


1 

0 
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€ 5 


Rtlritive perfor'''jrjiCG of ho’uri'-’tics of EDB considering 
numbtr of problens lor vdiich optiraobly nf the 
noftueuco is guaranteed (i,c. = 0 or 1 ) 


Problcn 


Ifeuriotic Ifo. 

1 2 3 


2x10 

2x25 


2x40 


3x10 


3x25 

3x40 


4x10 

4x25 

4x40 


5x10 

5x25 

5x40 


Total 
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Tatle 6 


Relative performance of various heuristics 
considering average number of tardy jobs 


Problem size 

Average No. of 

tardy 

gobs 

mxn 


Heuristic 

No. 



1 

2 

3 

4 

2x10 

1.8 

2.2 

2.4 

2.2 

2x25 

5.8 

5*6 

5.4 

5.0 

2x40 

12.4 

11.8 

11.8 

11.8 

5x10 

1.8 

1.6 

2.0 

1.6 

3x25 

5.0 

4.8 

4.8 

4.6 

3x40 

8.8 

9.4 

8.4 

9.6 

4x10 

1.2 

1.0 

1.0 

1.0 

4x25 

4.4 

4.4 

4.4 

3.8 

4x40 

5.8 

5.0 

5.4 

5.6 





5x10 

1.0 

0.8 

1.0 

0.8 


5.8 

3.2 

3.2 

2.8 

5x25 

5.6 

5.4 

5.8 

5.8 

5x40 




Overall Average 

No. of tardy jobs 

5.8 

4.6 

4.6 

4.5 
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Table 9 


Average 

execution times for 

problems 

of different 

sizes 

Problem 

]7o. of 

Average time 

token 

by 

size 

problems 

Heuris 

1 

Stic (in 
2 

milliseconds) 
3 4 

2x10 

5 

10 


10 


2x25 

5 

21 

18 

39 

19 

2x40 

5 

74 

90 

102 

72 

3x10 

5 

* 


14 

11 

3x25 

5 

32 

51 

34 

35 

3x40 

5 

32 

84 

93 

90 

4x:10 

5 

* 

14 


14 

4x25 

5 

41 

36 

49 

58 

4x40 

5 

78 

70 

72 

71 

5x10 

5 

12 

* 

13 

15 

5x25 

5 

40 

56 

60 

84 

5x40 

5 

83 

82 

121 

98 


* denotes the average execution time less than 10 milliseconds 
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Table 10 



Data for 

the 

ord ered 

flow 

shop 

problem Fo.l 


» 


1 

2 

Job k 

3 

4 

5 

6 





■^Ik 

5 

10 

12 

14 

15 

18 





■^2k 

26 

30 

55 

40 

50 

. 55 





^k 

150 

140 

150 200 

150 

300 








Table 

11 








Data for the 

ordered flov 

shop 

problem JJo , 2 





1 

2 

Job k 

3 4 

5 

6 

7 

8 

9 

10 

■^1k 

53 

32 

32 

31 

30 

25 

13 

12 

11 

10 


85 

75 

72 

70 

65 

60 

42 

40 

33 

33 

•^k 

150 

280 

390 

500 

250 

200 

220 

160 

140 

320 
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Table 12 

Data for the ordered flow shop problem Fo. 3 


M MM* wm 

1 

^2 

_ ^ 3 _ 

Job k 
4^ _ 

- - 5 _ 

- . 

7 


40 

42 

44 

46 

48 

49 

50 

■^2k 

50 

60 

70 

SO 

90 

95 

100 


160 

350 

300 

250 

180 

150 

200 
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01800 
01850 
01875 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
_ 02600 
102700 
,02800 
02900 
03000 
03)00 
032v)0 
'03300 
03400 
03500 
03550 
03600 
03700 
103800 
:03900 
04000 
' 04400 
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c 

f* 

c 


PRUflHA-! 1') 'il ,H(, fif 

r\i-i'.Y jn.ui K!M (--(vr ^-jub plcu' shop, 
t n'Mf ft (.••■rKGhPf A-2) 

:)rrK-:i^ V.h. M l.(Hf-) ,'fA\'(5) ,H(5) ,.'<TC5) 
r 1 '!!■ , 6, LI., T to ! 10,0! 50) 

l,pr-(bu),C(t.:b,L!50 

2) , i;iTAf.T(5Ct) ,1H50} , Jul'T, IT ( 50 ) , IJ.TAR ( 0 : 50 ) 

OA'l A I 'jPUT, luHT, lMiTl/2t!,2J ,22/ 
fiPK> (I‘9ITs;Ir.pHT,F.il,i.,=:'f-K,DAT' J 

nPKt.l OWIT=inUTl ,FILK='TT.OnT' ,RCCFSS='APPPMO' ) 

I T2. (JOT*, Access=' APPEND') 


fcir the explatioh of the \IARIABLE MAMES used PL. 


REF ER 


section 2.2. 


10 

r . 


28 

C 

c 

C32J 

CiU 


KEAU(JfsiPUT,4)AUOEPR 
OJ 331 EFFsl , JUuFPR 
HEAUdMPUT, 

on 10 j-i,H 

KEAtHlNPHT,=t=) (Til ,J) , 1 = 1 ,m) ,DtJ) 

PPCJ)=4 

CQ*iTI6UE 

SUBRUOTIfJE \;SRrpM ARRANGES THE 
i . EDO ORDER, 

CALL (0,H,PP) 

DlJ 28 R=i,N 
eDO(R) = PPCP.) 

NHITEC l(U(T,313) (PP(0) ,J=1,B) 

WRITE ( TOUT' , 32 3 ) C < I , t ( X , J ) , J=l , N ) # 1 = 1 1 (4 ) 
r.’P^ VC-.i ISa, LMB-NO. V20X, 50(12, 5X) ) 

5* i'-.' A‘‘ t 1 I ‘ , > 


•( 1 ) 


i . 

f' •>. 


(• J 'PL 
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C 

p 

c 

c 


(j-H.S'jirrl ! 


.5) 


’.r 

r . 

tnl'l 

E 1 L ( 5 ^ ' i'- 

I. i, r, 1,)0 

•j ! A<’= •'■r 
EOU- '«'= ' 

C,fVE ’’-lE I EU<' 

full ''ip'au-- EIL.'A T-.' '.<ET. SEQ, 
,'V APFi.t JI’H; -(■'!<!. K'S Gf,;- T !'(?!•' !, 

■.J<1 C''iG',= 1 ,4 

CALI f'T 5 o'C!' 

(FtrjLE.UT.I 14114 

1).1 2 1 ,s=l,f-; 

B I 



04500 
04600 
04700 
04900 
05000 
05100 
05200 
05300 
05400 
0550(1 
iibiiOn 
0 57ff0 
05800 
0590 0 
061 (H, 
0620U 

063 0 0 
07000 
07100 
07200 
07300 
07400 
0750 0 
07600 
07700 
07800 

07900 

07901 

07902 

07903 

07904 

07905 

07906 
08000 
08100 
08200 
08300 
08400 
08500 
08600 
08700 
08800 
08900 
09100 
09200 
09300 
09400 
09500 
09600 
09700 
09800 
09900 
10000 
10100 
10200 
10300 
10400 
10500 
10600 




20 


4C5f} 

JO'.U 

I 0 


I AS 

M 

i 1 1 
i'Jut 

4fol}0 
1 SO 


c 

c 

c 

c 

1 3U 

1 40 
165 


170 


625 

1010 

1030 


7 07 0 
1700 


1020 

C 

C1234 

40 

300 


PPCSJsKDlHO) 

6TAH=:Fni-;f.'j 
CALL CU'iPLi- 

w + 1 

'’.‘irf.ifv, -I 

If (-■'‘L.Lt;.! J 360,4 U{. 

CI'bL Clh\hLu 
!j'1 40 1 = 1, i>«;, 

iK ( (fc .,pp(3))-i irM(nn.or.o)io,4o 

‘i4X = 0 

i j f» ” 1 

IF (I.LP.l) 145,*i20 

f j < = 1 

G.> !•■! 707u 

IFt J ,h;;.2)3n4t5C» 

Lf; = i 

■3 J T't lr>5 

) 1’i( 4Mp'.>, i7 0i.,4bOo,4 90(‘i ICllDK 
CAf.f, n -If 2 
G1 16(1 
CALI. TT'^’CI 
Gn TLi 160 
CAta. rT2 
GG TO IbO 
CALL TOTAL 
DO 140 IJK=2,1 

IF Cn CPP(IJKl),GT,r.AX)130,140 


iUXsITCPPdJK) 1 
L K — I J K 
CUf»(TINUE 
PPtN+1 )=PP(LK) 

LCOMP=C(M,PP( ‘t + l 3 ) 

DO 170 A=LK,W 
PP(A3=PPCA+1) 

CALL COMPLH 
CALL TARDY 

IF C«T-LTAR)625, 625, 1010 

GO TO 1020 

DO 1030 BBB=t'i,LK+l,-l 
PPCB8B)=PPCHRB-n 
PP( +J) 

l^r (jg, 1^ 

PP€(Wt) = PPCLf<3 
00 1700 AsLK,^ 

Pis’ll'! -Vfi j 

CALu rij 'J-f, . 


6 1 . ; . 


M * 1 ' If 


6JT rS(H,i2i4 ) (PP(J) ,j3l,N3 
OX, 5(HT2, 

G t Tf) 2U 
C J.'. C f. aiK 
C'AiiL CO ^'FL ' 
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Paqe 



10700 

1 0 « 0 a 
10900 
11000 
11100 
U200 
lUOO 
11350 
1140 0 

1 1 s5 0-:} 

11 bOO 


I! 70>> 
UtfUU 
1 1 VOO 
12UHU 
1709O 
UlOO 
1220 if 

12iOu 
1240 0 
1250 0 
1270U 
I2yoo 
12900 
1300 0 
1305D 
13100 
13300 
13400 
1350 0 
13700 
1380 0 
13900 
14900 
15000 
15100 
15200 
15600 
15700 


15800 


C WRiTfCIOUT, 7893 (f ( v, ) , jj = i ,;.j) 

C789 PJk.-'A'I C10X,5O( 13, '-*)) 

Cl 2 f l^^.-ATCinx, IMi, , i J 1 1= 1‘= ' , 1 5 3 

Q\hh TARO 7 
GOTAR(KKK)= vf 

C .^RirR Cinijf-,2*>3KKiv,(FPCn,T = l,?0 

C25 8 3R ^fliTC IQX, M«'Ot-SKO, AT THK OF 

C I ribRAT-CO'' TS A.S Fn 

C 1 f.fi’j'^S'/lOX, 50(14)//) 

AKtiF c ii'iin ,55);'. r 
55 f JR ' A 1 (.1 UX, * IT = ' , 14) 

IF C '-7 35 10,55 0 

550 r-' Ah'CKFh-l ) )600,700,8y0 

-V > 1 ! F( l')‘i r, 1 so ) 

150 J,- S I ( 1 4,' < , * I' I^FFOljR , ,‘<. 1 , riF rAHl)7 

U-'ir S '-.‘'AS I Oi.-jA.Si' r>.& IT i.5 
1 ^ , r 4 } 

■; ) 1 ) 111 
C 70 i,f C i ')”(■ { oji-. 

70'' iij ru in 

sou OTAOsoT 

5l‘,J CO'lTlvUF 

111 CAOO « riME2) 

c , wRi rF(iouf,520){ppcn,i=i,*‘i) 

Cs2o F3RMAT(15X, MPT. KbO, IS AS FOLLOWS 

ICVIHX, 50(14)//) 

TIMF.3s:(HMF2.-TI^Fl) 

r£( lOUTl ,38)CUfrE,TlMF3 
88 F3R^.AT(10X, h,10) 

222 FOn.’^ATClOX, 'FUD OF pRo. i;0, — ',15) 

R9 CQNTINUF 

^RITECIOUl ,222)Fl'F 
WRITECIOUTl ,222)FPF 
3 33 Ca-TTlNOe 

500 8RITFCIDU1,595)' 

4.RITSCrUUTl ,5''3)' < •• FPn 

595 FORMATCiOX, *Eu.; ur ',15,' PRO, SESSIUi'l') 

39 CONTINUE 

SCOP 
Ell) 
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00010 

00011 


ApPfcJOlA B 

00012 

OUOl J 

\m 


00014 

r«| 

V# 


0001 b 


Paqe; 

0001 t,. 

c 

7 



Oi,iU1 h 

r* 


00 1 Hl;l 


rt.iiuns’ r 1 .'t 

0'?2‘'u 


f ’-n.icr r t -i <.- 2 ) 



c ) . i'!ii?ovo,TMh .F:s,n,>«,i.i.,TC-5;io,o;5u) , 



IPPtbvj ,Cii-'^i(o;b,<):bo 

00 40 0 


2 J , M ! Al. f CbO ) , A( S'.' 1 - lUOT 

OOHOO 


O 1 S KA = 1 , '-1 

0 0 0 0 U 

s 

C I 1 ) J=CU..ptKK-l,PPCU)+TCKB ,PPC1 ) ) 

00 700 


OOBOn 


•>| t:i :<si,4 

OO'hOU 


o.i 10 J = 2 ,a 

0 i 0 0 U 


CK£C '1 V(K,PP(J-1 ) ) .Gf. cn>;p(K;-l ,PPCJ) 3 3 300,400 

0 1 i 0 0 

J'MJ 

C }'’PtK.PPCJ))sCti^iP(K,PPiJ-l ) J+TCK,P?(J) 3 

G 1 n 1 0 

01 20 0 


01300 

100 

C ) •ipCK,pptJ))=:CnMp(K-l,pPCJ) )+nK,PPC.J)3 

01400 

10 

cn ’'jrt'vue 

01500 

fy 

wHrre:£iuuiM2 3Cnpp(M,pp(f;2) 

FOPf'A'i (lOX, 'CU«PL\,-TIfc!K, 1 i“'*=',I53 

Olt.OO 

zv?. 

01700 


PfiTllPO 

01 H0() 


E 4D 

0190v) 

IW> AiC m M 


02000 

*4, m m m m 


02100 


SUBPOOTI^e TAPO-7 

0220U 


T, 3 P H C T. T 1 r.- T E r, E ( A - Z 3 

02340 


'> i .'PhES ,« , LL , r ( -5 no , 0 : 50 ) , 

lPPC50),r.;' .. Cu;f,0;50 

13 ,TQT-: r(s: ) 4503 ,iniiT 

023B0 


02400 


02500 


TAfti'.Y''' = 0 

02000 


^:o=•' 

02700 


■! 1 1-' 1 = 1, W 

02800 


IF (tCUMPtM.pPCl))-.U(PP(I J) ) .LE. 0)10,20 

TAROYwsTAROY.N + l 

02900 

20 

03000 


TAR7ESsTARNES+CC0MPCH,PPCI) )-DCPP(I333 

03100 

! 

c . 1 n 

03200 


^ -l) 1^(4 

03300 


P.IP ..iUl ox, MO, OF 'lAROY J06SV16X,I43 

03400 


i.Ui;r,b':OT4RKES . „ . 

03500 


41 (?(.•>:, '-“I'A'rAT. TAPDV' ESS FOR ABOVE SEQ, IS a*,l4/) 

03600 


■'XITEC lii )'j’,hb5) Cl. If PP("')) 

03700 

*3 ''3 ^ h 

!^3.< Al{!5/.,* C'i-.Fl.i- ri.' ^ TIME s', 133 

03800 



03900 


. : ) 

04000 

^ W M M 

4, 


04100 

^ w Ifw 


04200 


SUHp.UUTI.iE irs 

04300 


J-lPLICIt IwTEGfcPCA-2) 

04340 


CO f AHDY t , Trtpf. ES , M , 8 , LI, , T C -5 n 0 , 0 : 50 3 , 

04380 


IPPtSO) ,C0>;F£0:6,usb0 
l),10TALt(50),r>(50),lOUT 

Bl i 

04500 




04600 
04700 
04800 
04900 
05000 
05100 
U5200 
0b30tJ 
( ) S ^ o) 


un 10 j=i,2 
T.jrA0TCP8{J)) = 0 
Da 10 Ksi.M 

10 TUTAf.T(PPCJ) )=Tf'‘rAL*r(PP(J) )+Jr({«.,FP(J) ) 

C ir)UT,25) (Tr-TAI,TCPP(J) ) ,J=1 ,'0 

C25 FJiOviA-i (4 OX, 'TUTaI-TIA'P'/IUX ,50( 14, ) 

t. 0 
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ftPPRtJDlX P 


OOlOO 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

OtlOO 

01200 

01300 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 


SUBROUTINE TTBCl 
IMPT 4 TCTT IMTEGERCA-Z) 

COMMON TAROYN-TARNES. M,N,LO,TC-5UO 
tPPC503,CUMPC6s6,0s50 
l),TnTALT(50),DC50),IQUT,ITC50l 
DO 10 0=1, N 

10 ITCPP(J))=TC1,PP(J)) 

RETURN 

END 

mf m m m mm m m m 


SUBROUTINE TTMC2 

IMPl^TCTT INTEGERCA-Z) 

common TARDyN,TARNES,M,N,H.,TC-5n0 

1PP(50) ,C0MPC6:6,0S50 

1) ,TOTAL,TC50 3,DC50) ,IOUT,ITC50) 

DO 10 J= 1 ,N 

10 1TCPPCJ))=TCM,PPCJ)) 

RETURN 

END 
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SUBROUTINE IT2 

IMPLICIT INTEGERCA-Z5 

COMMON TARDYN,TARNES£M,N,LL,TC-5S10 

IPPC50) ,COMPC0l6£0!5O 

l),TQTALTC50),DC50),IOUT,ITC50) 

DO 10 J=1,N 

10 ITCPPCJ))sCCM,PPCU))»C(M,PP(J-1)) 

RETURN 
END 
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c 


mwWmMxr'fmWMMayfffivOTMwafWM 


0J501 , 


OlSOl, 


0*50), 


B21 



APPji'-hJT C 


aolou 


p^nr;, r) £;k ifra-if: PArorif" ijat.\ por,:-,/*.! n-c case. 

00200 


r tPLlClT fA-Z) 

00300 


h(5uj0) 

00400 


nr TOTATi( I'lOO) ,T( 5, tuoo) ,0(10001 

00450 


•OM’A TuOI ,1 vPUJ , £..C9,10 mTI/21 ,20,0,22/ 

0070(* 


OOS- (?l .1 ( = 1 'Ip'ir,! n.Ks'Hl .DA A*) 

00H0-; 


'IPK' ( .r- 1 islf’.ij i f,e,r»? \ 1 .Hm f ') 

00i-i5 ! 

‘’I 

4 .* 

( < ! '.sr.'Mj'i ,h ti.Fs;' i i -.F.our^ACCKftSs'APt’FwO' ) 

AOOOu 

m 

TAf^T, A ri .F.( 1 l-.KI ) 

COOP'S 


p-; ‘i'.jf 1 ir'i'i', 9) 

00‘35'.‘ 


-1! 9 1 11 XV'Z=1 SPSn.j 

01000 


R-:Ai (1 ,*) ■otrti.p.. ,-,r,]ZKl ,MSTZK2,TSIZF,1,TSIZF:2 

1)1 lOt, 


1 , LC-ii’i.p, uj ri;?z,oi:,T IV 

01 bOO 


1 r )i p ?=( 1 1( isiZi-2-TsrzFi 3/i3Ti-:P)» 

01t>40 


1(U( ST ■^K2-^;6lZF1 I)f:r-:1FP« 

01o50 


'RHlrU M'll 'StSuM 

01700 


-iKr»Fc rnoT, ixoTpp 

01 BOu 


IVJ 200 -s.lSlZKl ,VSTZK2 

01 90v 


m ISO :v=tsrzFi ,i6rzF:2,i5fFP 

02000 


CtKUi, hA.>0n»'>CR,R3 A’f'FX, INCF ) 

02100 


0.1 lou jjj=i .nonFFR 

02200 


-vHnf-:ciu'jT,20i3 9,''; 

02i00 

Pul 

F JK/CATCIOX, I4,5X,I4J 

02400 


OD 20 0=1,/ 

02500 


0-J 20 X = l,\ 

02600 

1 1 1 

Pl.«f-FX = HTf..OFX + l 

02700 


ifchimikx-=5du 0) CAOt r; ' o/n-»iCri,KX^n)i:x, uiCH) 

02B0xj 


It<=:RCPT.''’0KX)4t00 

02900 


IF (IF. 1 - 0 . 0 ) yn Tu in 

03000 


rFCiH.yT.LTvi'iM) bo TO ni 

031O0 


It J, U = 1H 

0320U 

20 

C'-J^TT^OF, 

03300 


oa 3U 1=1, N 

03400 


T3TAL( I)=il 

03500 


on 22 J=1,M 

03600 

22 

TnTAr,i(I) = TOTA15Cl)V!'( J,T) 

03700 

25 

RIMf)SX=KlMDEX+1 

0380U 


lFCRTFDfc;X ==SOOO) CALI. RAODn!4(H,ftlNDEX, INCH) 

03900 


IRaH(m«!)EX)1'5000u 

04000 


tF C.£R==0,nR.l9.(;v,i,lMDUK.OK,lR*L'UTOFA|,CX) ) GO TO 25 

04100 


Dtr.) = Ift 

O42O0 


WRitFcinae,*) (Tco, r),j=i j 

04300 

30 

C'i Hi'iuit; 

04400 

100 

C ) UMlZ 

04500 

150 

C: 11 Milt; 

04600 

200 

C7 U-Uti, 

04650 


c . -T-rrrT '■' 2 ) 

04675 


’f* * - . J ' 2 1’ ! 1 c 1 

04687 

c 

lOUTl ,993TIMEi 

04690 

9999 

C IV M'-AF 

04693 

99 

T !< '-X , 'TIME-TAKEN =',16) 

04700 


3 - np 

04800 


F 'j D 

04900 



05000 

C 


05100 


S JO-nnrf F 04 .CP,R1FDBX,IMCB) 

05200 


,)t i 4fSJ001 

05300 


r . f-FO./i n X, 1-CR,15ED 

^ ^ A 



i)5b0a 
UbbOO 
Ob 700 

05 BOO 
0590 n 
060 00 

06 J 00 


! 'Ci'=J OCI- fi 
r.OKOsi5^;')i-i2? + i."CB 
jsSO )M 

CAt,{. t.’G00(TSKin'.,7O 
01 jn;>:!(=i 
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